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Abstract 

We present a method for finding lines of constant physics in the confinement phase 
of the SU(2) Higgs model on the lattice. The model is considered at finite values of the 
cut-off where it behaves like an effective field theory with three independent couplings. 
In particular, a renormalised quantity sensitive to a variation of the bare Higgs quartic 
self-coupling is constructed from generalised Binder cumulants. Numerical results for 
the non-perturbative matching of the bare parameters of the model between (3 = 2.2 
and (3 = 2.4 are presented. 
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1 Introduction 

The conventional bare parameter space of the SU(2) Higgs model jlj is defined by three 
parameters: = 4/g 2 , where g is the bare gauge coupling, the hopping parameter k 
and the bare Higgs quartic self-coupling A. The continuum limit of the lattice model 
is at (3 = oo and k = 1/8. In general, the physics in the continuum limit of a lattice 
regularised field theory with several relevant couplings depends on the way in which 
this limit is approached. The approach to the continuum is represented in the bare 
parameter space by lines on which dimensionless and renormalised ratios F of physical 
quantities are constant. These lines are called lines of constant physics (LCP) Q. 

It is well accepted - supported by the weak gauge coupling expansion and by 
numerical simulations - that the scalar part of the SU(2) Higgs model in four space- 
time dimensions is a trivial theory Q , which means that in the continuum limit of the 
SU(2) Higgs model the renormalised scalar coupling Ar is always zero and therefore there 
are two relevant couplings only. Nevertheless, since the coupling A is only marginally 
irrelevant, the SU(2) Higgs model can be considered as an effective field theory with 
three non-zero renormalised couplings at finite values of the cut-off (the inverse lattice 
spacing a). In this article, we describe a method for studying the scaling properties of 
this effective field theory. 

In practice, the LCPs in the SU(2) Higgs model are constructed by changing (5 
(which corresponds to a change in the lattice spacing) while renormalising the parame- 
ters n and A in order to keep two independent physical quantities F\ and F2 constant: 

Fj(/3,K,A) = constant, £ = 1,2. (1) 

Along a LCP, the value of a general physical quantity F (different than F\ and F2) 
shows a dependency on the lattice spacing like ||]J| 

F(a)\ FuF2 = F + o((a/r ) 2 ) , (2) 

where F is independent of a and ro is the physical scale, like the Sommer scale ||. 
The O ((a/ro) 2 ) terms arise from irrelevant operators and are called lattice artifacts. 
The absence of O(a/ro) corrections on the right-hand side of eq. (Q) for a scalar theory 
is well accepted. Because the renormalised coupling Ar vanishes only logarithmically 
along a LCP @], i.e. like l/ln(ro/a), it is possible to find a range of values of the lattice 
spacing where the corrections on the right-hand side of eq. (||) are negligible but Ar is 
still sizeable. In this range, which we refer to as the scaling region, for a given lattice 
spacing the bare coupling A can be tuned such that Ar assumes a particular value: the 
effective field theory described by the SU(2) Higgs model with finite cut-off has three 
independent couplings. 

The motivation behind the work presented in this article is the following. In a 
gauge theory with matter fields in the fundamental representation of the gauge group, 
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the potential between a static quark and a static anti-quarl^] is expected to flatten 
at large separation of the static sources. This phenomenon is called string breaking 
||. Whereas string breaking has not been observed in QCD simulations yet pUlC|l, it 



has been established by numerical simulations of the four-dimensional [11] and three- 
dimensional (T^j SU(2) Higgs model, where a new method for determining the static 
potential has been used. The phase diagram of the SU(2) Higgs model || is divided 
into two "phases": the Higgs "phase" at large values of k and the confinement "phase" ^ 
at smaller k. In the confinement "phase", the model has confinement properties like 
QCD 0,13] and is therefore suitable for testing the method of determining the static 
potential. In order to investigate the size of scaling violations in the measurements of 
the static potential a way to define LCPs in the confinement "phase" of the SU(2) Higgs 
model has to be found. We anticipate that our results for the static potential in the 
four-dimensional SU(2) Higgs model [|l4|,|ll| show compatibility with scaling within tiny 
errors already at moderate (3 values. 

In three space-time dimensions, the LCPs can be exactly determined within pertur- 
bation theory [16] due to super-renormalisability of the theory. In the four-dimensional 
model, which we consider in this article, the variations of k and A along the LCP have 
to be determined non-perturbatively. We construct two physical quantities, F\ and F2, 
which have a sufficiently independent sensitivity on k and A in order to be taken for 
the definition of a LCP. The quantity F\ is constructed from the static potential and is 
sensitive to the mass of the Higgs field which determines the string breaking distance 
?"b- The quantity F% is a generalised Binder cumulant constructed from the connected 
p-point functions of a zero-momentum gauge-invariant Higgs field: it has proven to be 
sensitive enough to the Higgs quartic self-coupling that we can use it to construct a 
LCP. We present numerical results for the renormalisation (or matching) of the bare 
parameters along a LCP between (3 = 2.2 and = 2.4. 



2 Definitions 



A sensible choice for the physical quantities F\ and F2, which allows us to match the 
parameters k and A along a LCP, would be to find an F± which is strongly dependent 
on k and an F2 which is strongly dependent on A. In this sense, a good quantity Fx can 
be defined as 



Fx = r [2fi - Vb(r )] , (3) 

where Vo(r) is the static potential and \x is the mass of a static- light meson, the lowest 
bound state of a static quark and the light dynamical Higgs field. The difference 1\x — 
Vo(r) is a renormalised quantity free of self-energy contributions of the static sources, 

1 A static quark (anti-quark) is an external source in the (complex conjugate of the) fundamental 
representation of the gauge group. 

2 At small values of (3 there is an analytic connection between the two "phases": therefore, we put 
"phase" in quotation marks. 
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which diverge in the continuum like a -1 . Because the static potential flattens reaching 
its asymptotic value lim r _ >0O Vo (r) = 2/x, the quantity F% eq. (||) is strictly related to 
the string breaking distance r^, around which the static potential starts flattening out, 
and is mainly sensitive to the value of fi and hence to k, which determines the mass of 
the Higgs field. 

A remark should be added concerning finite size effects. Because physics shows a 
dependence on the physical size of the system, the latter is part of the definition of a 
physical quantity F. Along a LCP, the physical size of the lattice L/vq has either to be 
kept constant or at least large enough to avoid finite size effects. As concerns the static 
potentials and the static-light meson spectrum, the relevant quantities to compare are 
the lattice size L and the mass gap in the gauge invariant sector of the theory, the Higgs 
mass rajj. If Lmn 3> 1 finite size effects are expected to be negligible. This condition 
is fulfilled in our numerical computations |jl4j| . 

Now, we are looking for a physical quantity F2 which is sensitive to a variation 
of the parameter A, once the parameter k is matched by using iq eq. (^). There exist 
perturbative equations |7j describing the change of A and g along a LCP in the weak 
coupling limit. But the use of these equations is not justified in our case because the bare 
couplings are not small enough. The change of A along a LCP has to be computed in 
a non-perturbative way. Exploratory results |ij indicated that the physics of the SU(2) 
Higgs model in the confinement "phase" is weakly dependent on A once a physical 
quantity such as F\ eq. @ is kept fixed. Our studies of correlation functions from 
which the Higgs and W-boson mass can be extracted confirm these observations [14]. 

Here, we consider a zero-momentum gauge-invariant Higgs field 

s = ^£ 5 (*) = -E^c^c*), (4) 

where &(x) is the complex SU(2) doublet Higgs field and f2 is the number of lattice 
points. The field S(x) is a composite field. The renormalisation of composite fields is a 
complicated issue and we refer to textbooks, e.g [pl|Jl^1 , for a detailed discussion. The 
field S(x) gets renormalised like Sr(x) = Z\ + Z$S(x). If one considers the connected 
p-point functions G Cj r(xi, x p ) = (Sr(xi) • • • 5 , r(x p )) c , the additive renormalisation 
Z\ cancels and the multiplicative renormalisation Zs is chosen such that the functions 
G c ,r(xi, x p ) have a well-defined continuum limit, provided the arguments x±,...,x p 
are kept at non-zero distance from one another [fl9| . 

Additional care is required if some of the arguments Gc,r{xi, Xp) 

coincide, because in this case divergences can arise in taking the continuum limit. This 
situation can be analysed in the continuum with the help of the operator-product ex- 
pansion of Wilson [pOl , according to which 



(SrOi) • • • S R (x k )) c ^ ~ |^2fc (e -> 0) , (5) 

in the limit of all coincident arguments. The power of the divergence on the right-hand 
side of eq. (J5f) corresponds to the naive dimensional counting. When we take the fields 
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Figure l: The cumulants C3 and C4 defined in eq. (||) as functions of the lattice size at 
(3 = 2.4. 



Sr(x) in eq. (pi) at zero momentum we get the expression 
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vp- 



d 4 xi ■ • • d 4 x p (5 R (xi) • • • S R (x p ))c 



d 4 yi • • • d%_i (Sr. (yi) ■ • • S R (yp_i)S R (0)) c 



(6) 



where V is the physical volume of the torus and in the second line we used the translation 
invariance property of G cR {x\, x p ). From eq. (||) and naive dimensional counting, 
we infer that in the limit when yi — » (i = 1, ...,p — 1) in eq. (^) the integral is finite 
except in the case p = 2 where it is naively logarithmic divergent. This problem can be 
cured by defining a modified connected 2-point function 



(4 2) } c = i J d% (S R (y)S R (0)} 



, 2/0 
c sin I — 7T 



(7) 



where T is the time extension of the torus: we use T = L = V 1 / 4 . The sine function 
gives a contribution proportional to yo when y — > and therefore it makes the integral 
finite. In order to be sure that all divergences associated with coincident arguments 
in eq. (||) are regularised for p > 2, it remains to consider the case of two coincident 
arguments, e.g. y\ — * yi. In contrast with the case p = 2, there are still two integrations 
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Ak 
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K 


Ak 


0.4 


0.2615 


0.0001 


0.757 


0.2973 


0.0001 


0.5 


0.2737 


0.0001 


0.96 


0.3108 


0.0001 


0.55 


0.2792 


0.0002 


1.0 


0.3131 


0.0001 


0.7 


0.2929 


0.0001 


1.5 


0.3351 


0.0002 



Table l: Values of k for different A's determined at (3 = 2.2 by using the matching 
condition eq. dTT|). The values of Ak are the uncertainties associated with the matching. 



over yi and y2 and the integral is finite. We can now define generalised Binder cumulants 

CP = 3, 4,...). (8) 
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From the above discussion, it follows that the couplings c p are well-defined and renor- 
malised. Notice that the bare fields can be used in the definition of eq. (||) because the 
multiplicative renormalisation factors Z$ cancel. 

The cumulants c p can be defined on the lattice using the same expression eq. (||): 
the field s is given in eq. @) and 

(s {2) ) c = ^E(S(x)S(y)) cS in(^^^ (9) 

is the lattice version of eq. (|7j). In the next section, we present Monte Carlo results of 
the computation of the cumulants c p on the lattice which show a linear sensitivity to 
a variation of the parameter A, once the physical quantity F± eq. (|3|) is kept fixed. We 
define the second physical quantity F 2 characterising the LCP to be 

F2 = c 3 | L/ro , (10) 

where the physical size L/ro of the lattice on which C3 is computed has to be kept con- 
stant along the LCP. This is required because the cumulants c p are finite size couplings, 
as we are going to see in the next section. 

Our final result is that physics in the confinement "phase" of the SU(2) Higgs 
Model can be reproduced at different values of the lattice spacing by changing the bare 
parameters k and A such that the physical quantities F\ eq. @ and F2 eq. ([ji]) are kept 
constant. In the next section, we construct two points on a LCP at (3 values 2.2 and 
2.4 by computing Fi and F2 in Monte Carlo simulations. 

3 Numerical results 

At /3 = 2.4, we take the parameter set k = 0.2759, A = 0.7 Jl4| and we determine the 
corresponding values Ff and F£ (referred as matching conditions) of the physical quan- 
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Figure 2: The curve k(A) at /3 = 2.2 from the matching condition eq. (11). The solid line 
represents the polynomial fit eq. (|l5|) . The vertical dotted band represents the matched 
value of A (with error) along the LCP eq. (|l7|) . 



tities F\ eq. (||) and F2 eq. (|i~0|) that define the LCP. From the quantity Fi, computed 
on a 32 4 lattice, we obtain the first matching condition 

Ft = Ff = 1.26, (11) 

with a statistical error AF* = 0.02. 

The cumulants c p defined in eq. @ are very small for large values of p: in the 
Monte Carlo simulations, we were able to obtain a significant signal up to p = 5. In a 
massive theory, from the translation invariance property of connected p-point functions 
(see eq. (0)) the cumulants are expected to scale for large volumina V of the torus like 

/ , \ l-p/2 

c p oc (tt&V) , (12) 

where mn is the Higgs mass. We computed C3 and C4 at /3 = 2.4 for different lattice 
sizes: the results are shown in Fig. |l|. We observe strong variations and therefore, a 
precise matching of the physical lattice size L/r$ along the LCP is needed in order to use 
F2 = C3 for defining the LCP. We choose to compute C3 on a 7 4 lattice at /3 = 2.2. The 
physical lattice size is small enough to obtain a good signal and the matched lattice size 
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Figure 3: The A-dependence of the cumulant C3 at (3 = 2.2. The values for the parameter 
k are the same as in Table EL The horizontal dotted band represents the second matching 
condition eq. (|l3|). From the slope of 03(A), we determine the matched value A*. 

at (3 = 2.4 is L/a = 13.3(3)0 in a region where C3 has no strong variations (see Fig. EJ). 
Interpolating the values of C3 at /3 = 2.4 shown in Fig. EJ to the matched lattice size 
gives the second matching condition 

F 2 = F 2 * = 3.17, (13) 

with an error Ai 7 ^ = 0.07. 

Now, we are ready to match the parameters k and A at j3 = 2.2. We start with 
matching k for different values of A by computing F\ on a 20 4 lattice and requiring 
eq. ([n]) to be fulfilled. The resulting (A, k) pairs are listed in Table [l| and plotted in 
Fig. ^. We fit these points with a polynomial 

M 

k(X) = ^a fc (A-l) fc . (14) 

k=0 

A good description of our data is obtained for M = 4: 

k(X) = 0.3131 + 0.0564 (A - 1) - 0.0286 (A - l) 2 

+ 0.0198 (A - l) 3 -0.0246 (A - l) 4 , (15) 
3 The error comes from the determination of rg/a at the two /3 values. 
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Figure 4: The A-dependence of the cumulant C4 at (3 = 2.2. The values for the parameter 
k are the same as in Table |]. The horizontal dotted line is the value of C4 at j3 = 2.4 
for only approximately matched physical lattice size. The vertical dotted line is the 
matched value A* along the LCP eq. (|l7|). 



and a comparison of this expression with the data points is shown in Fig. 0. The 
deviation between the fitted curve eq. (jl5|) and the data points is smaller than the 
statistical errors of the data. When k is evaluated using eq. (]T5|), an uncertainty of 
1 • 10~ 4 should be assigned in the range 0.4 < A < 1.06 growing up to 4 ■ 10 -4 in the 
range 1.06 < A < 1.5. As a by-product of the simulations for the matching of n we 
obtain for the scale ro the typical value r^/a ~ 2.8 when F\ takes values near F*. This 
has to be compared with the value r^/a = 5.29(6) obtained at (5 = 2.4: the resulting 
change in the lattice spacing is 



a{(3 = 2.2) 



a((3 = 2.4) 



= 1.90(4) , (16) 

Fi=FT 



where at (3 = 2.2 we take A = 1.0. 

The matching of A at = 2.2 is performed by computing F2 = C3 on a 7 4 lattice 
taking the parameter points (A, k) of Table ||. The results are shown in Fig. [|. The 
matched value A* is determined by requiring the matching condition eq. ji~3| ) (repre- 
sented in Fig. ||]by the horizontal dotted band) to be fulfilled: we obtain A* = 0.96(10), 
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where the error is determined from the approximate slope of 03(A). Summarising our 
numerical results, the points 



(3 = 2.2, k* = k(X*) 7 A* = 0.96(10) 

and (17) 
= 2.4, k = 0.2759, A = 0.7 



lie on the LCP defined by the matching conditions eq. (11) and eq. (|13|). The value of 
K at /3 = 2.2 is determined using the fit eq. fll5|). 

In Fig. |], we show the behavior of the cumulant C4 on the LCP eq. (|l7]). The 
horizontal dotted band is the approximate value of C4 at /3 = 2.4: the physical lattice 
size is only approximately matched. The values of C4 at j3 = 2.2 when A is in the 
matched interval around A* are compatible with scaling. 



4 Conclusions 



We have described a non-perturbative method for matching the parameters k and A 
along a LCP in the confinement "phase" of the SU(2) Higgs model. We consider the 
model at finite values of the cut-off, where it describes an effective field theory with 
a non-zero renormalised Higgs quartic self-coupling. The lattice spacing is varied by 
changing (5 while two physical quantities F\ and F2 are kept constant. A require- 
ment for F± and F2 is that they have a sufficiently independent sensitivity to k and A. 
The first quantity F± is derived from the static potential and is mainly sensitive to k. 
More problematic is the definition of the second quantity F2 because the physics in the 
confinement "phase" is weakly dependent on A. We construct couplings, which are gen- 
eralised Binder cumulants, from the connected p-point functions of a zero-momentum 
gauge-invariant Higgs field. In particular, we define a modified 2-point function in or- 
der that no divergencies arise in taking the continuum limit. In the range of A that we 
investigated numerically, these couplings show a linear dependence on A that allows for 
the matching. We determine by Monte Carlo simulations a LCP at (5 values 2.2 and 
2.4. Based on these results, it is possible to study the scaling properties of physical 
quantities such as the static potential |i"5|| . 
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